Suppression of the high p T charged hadron Raa at the LHC 
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We present a parameter free postdiction of the high-pr charged-hadron nuclear modification 
factor (-Raa) in two centralities, measured by the CMS collaboration in Pb-Pb collisions at the 
Large Hadron Collider (LHC). The evolution of the bulk medium is modeled using viscous fluid 
dynamics, with parameters adjusted to describe the soft hadron yields and elliptic flow. Assuming 
the dominance of radiative energy loss, we compute the medium modification of the Raa using a 
perturbative QCD based formalism, the higher twist scheme. The transverse momentum diffusion 
coefficient q is assumed to scale with the entropy density and normalized by fitting the Raa in 
the most central Au-Au collisions at the Relativistic Heavy-Ion Collider (RHIC). This set up is 
validated in non-central Au-Au collisions at RHIC and then extrapolated to Pb-Pb collisions at 
the LHC, keeping the relation between q and entropy density unchanged. We obtain a satisfactory 
description of the CMS Raa over the pr range from 10-100 GeV. 
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Jet quenching, or the modification of hard jets in 
a dense extended medium, is one of the most stud- 
ied discoveries at the Relativistic Heavy-Ion Collider 
(RHIC) pQ. Numerous experiments have established 
the suppression of high transverse momentum (high px) 
hadrons [2] at RHIC energies. There now exist volumes of 
theoretical calculations based on both perturbative QCD 
(pQCD) [3] and the AdS/CFT conjecture 14] aiming to 
describe jet modification at RHIC. Collisions of heavy- 
ions at the LHC, as measured by all three detectors (AT- 
LAS, ALICE, CMS), have demonstrated significant evi- 
dence for the modification of hard jets [5H7]. 

Before any such measurement may be applied to the 
detailed study of the structure of the dense medium, the- 
oretical calculations must be able to describe the basic jet 
quenching observables. A standard jet quenching mea- 
sure, as established by the RHIC experiments, is the nu- 
clear modification factor Raa, defined as, 
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In the equation above, d 2 NAA{b m i n , b max ) represents the 
yield of hadrons in narrow bins of px, rapidity (y) and 
centrality (designated by a range of impact parameters 
bmin to b max ) of the heavy- ion collision. In the denom- 
inator, Nbm represents the number of binary nucleon- 
nucleon collisions in the same centrality bin and d 2 N pp 
represents the yield of hadrons in p-p collisions, in the 
same bin of pr and y. So far both ALICE and CMS have 
reported the Raa of charged hadrons at the LHC [7] . 

In this Letter, we present a parameter free comparison 
with the 0-5% and 10-30% CMS data 0. The calcula- 
tion consists of two factorized parts: A fluid dynamical 
simulation whose initial conditions and transport coeffi- 
cients were tuned to describe the hadron spectrum and 



elliptic flow at pr<2 GeV in Au-Au collisions at RHIC 
and successfully extrapolated to Pb-Pb collisions at the 
LHC [SJ, and a pQCD based radiative energy loss calcu- 
lation which computes the medium modified spectrum of 
high pt hadrons. Unlike soft hadrons, these hard hadrons 
will be assumed to stem from the fragmentation of hard 
jets, modified due to passage through the medium. 

By "parameter free" we mean that the transverse 
momentum diffusion coefficient q, which controls the 
amount of radiative energy loss encountered by a hard 
jet, will be related to intrinsic quantities in the fluid dy- 
namics simulation which are entirely controlled by the 
soft hadron obsrvables. In this case, it will be scaled 
with the temperature dependent entropy density s: 



q(r,t) = q s(T{?,t)) , 
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where qo refers to the value of the transport coefficient at 
the highest RHIC temperature of Tq = 344 MeV, reached 
at the center of the 0-5% central Au-Au collisions at 0.6 
fm/c, and v± refers to the local flow velocity transverse 
to the jet. Once related to temperature in this way, q is 
completely controlled by the fluid dynamics simulation 
at any point in space-time, in any collision at any energy 
or centrality. Given Eq. ([2]), the space-time dependent 
entropy density in the simulation for LHC energies will 
predict the value of q at any space-time point. 

The remaining Letter is organized as follows: after a 
brief review of the essential features of the higher twist 
method to describe parton energy loss, we discuss spe- 
cific ingredients of the current calculation. This will be 
followed by a comparison to Raa in central and semi- 
central collisions at RHIC (i.e. Au-Au at y/s = 200 
GeV/nucleon), as well as the in and out of plane Raa 
in semi-central RHIC events; go will be dialed to obtain 
a best fit to this data. Following this the Raa will be 
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computed for 0-5% and 10-30% centrality events at the 
LHC (i.e. Pb-Pb at y/s = 2.76 TeV). We will conclude 
with a discussion of possible reasons why the theory un- 
derpredicts the Raa at Pt < 8 GeV at the LHC. 

It will be assumed that processes such as the produc- 
tion of high px hadrons which engender a hard scale 
(Q 3> Aqc£>) through much of their space-time evolution 
may be computed using pQCD. Furthermore, the hard 
scale allows for the use of factorization, effectively sepa- 
rating soft subprocesses from the part where the scales 
involved are hard [5] ■ In such a factorized formalism, the 
invariant cross section to produce a high pr hadron in a 
heavy-ion collision may be expressed as 
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In the equation above, y,pr represents the rapidity and 
transverse momentum of the detected parton, K is a 
multiplicative factor to account for higher order cor- 
rections, r is the location of the jet production ver- 
tex, tAjB(r ± b/2) is the nuclear thickness function 
[tA(r + b/2) = J dzp{f +6/2, z) with p(f,z) the nu- 
cleon density inside a nucleus], x a (x b ) represents the mo- 
mentum fractions of the incoming partons, Q 2 ~ pin 
is the hard scale of the process, G A j b B (x a , Q 2 ) is the 
nuclear parton distribution function, da ab ^ c d/dt is the 
short distance cross section for incoming partons a, 6 to 
scatter and produce partons c, d with a squared momen- 
tum transfer t = (p c — p a ) 2 , D represents the medium 
modified fragmentation with hadronic momentum frac- 
tion z = pt/PcT' and (L(r,0j) is the distance travelled 
by a jet produced at r and propagating at an angle 9j 
with respect to the reaction plane. Both 6 and f are two 
dimensional vectors transverse to the beam direction. All 
calculations will be carried out at mid-rapidity (y = 0). 

In this effort we will ignore shadowing corrections, both 
at LHC and at RHIC. This will give the calculated RHIC 
Raa a slightly rising slope with pr, compared to previ- 
ous calculations. The medium modified fragmentation 
function is calculated by solving an in-medium DGLAP 
evolution equation [10] valid when Q 2 3> qrj (where, 
Tf ~ p/Q 2 is the lifetime of a parton with energy p and 
virtuality Q 2 , undergoing a split), 
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where, P(y) is the Altarelli-Parisi splitting function for a 
parton c to split into two other partons which carry y and 
1 — y fractions of the momentum of c. The in-medium 
single emission kernel KI' 8 q 2 [Tl] is given as 
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involving the space-time dependent transport coefficient 
q introduced in Eq. pi. 

As stated in Eq. ||2j, the jet transport parameter de- 
pends on the local entropy density s. The entropy density 
is obtained from a 2 + ID viscous fluid dynamical simula- 
tion where the input parameters, such as the magnitude 
of the components of the initial energy momentum tensor 
as well as the viscosity and the final freezeout criteria are 
tuned to obtain the best fit with the spectra and elliptic 
flow of hadrons with pr < 2 GeV. As this Letter will 
deal solely with the spectra and azimuthal anisotropy at 
high px (pt > 6 GeV), we will only provide the most 
salient features of these simulations in the current paper 
and direct the reader to Refs. [T2J H3j for further details. 

In the fluid dynamical simulations used, the ther- 
malization time is t = 0.6 fm/c for both RHIC and 
LHC collisions. Initial conditions are based on the KLN 
parametrization of the collision of two nuclei with satu- 
rated gluon distributions [12] ■ The jets are assumed to 
be produced according to a binary collision profile at the 
moment of collision, i.e. at t = 0. For times between 
to and To, the conditions at tq are extrapolated back to 
to = fm/c unchanged, i.e. we assume that the soft 
medium remains unchanged from to 0.6 fm/c. 

Like prior calculations of jet modification in a fluid 
dynamical set up, the hadronic phase is also simulated 
as a viscous fluid and not in terms of a hadronic cas- 
cade. While such hybrid calculations are currently avail- 
able, we have picked the simplest simulation for this 
attempt to calculate the Raa at the LHC. As we use 
the entropy density to scale the local value of q, there 
is no extra rescaling factor for the hadronic phase as 
in prior attempts where q was scaled with T 3 [14j [15] . 
The produced jets are assumed to decouple from the 
medium when the local temperature reaches 160 MeV. 
Given the above conditions, a space-time profile for q is 
obtained for jets traveling in all directions, starting at 
t = at any location in the medium, and vanishing when 
T(r, t) = 160 MeV. This profile is completely specified 
with the value of go = 2.2GcV 2 /fin at T= 344 MeV which 
is the highest temperature reached at RHIC, at f = 
in the 0-5% most central collisions at t = To. This qo is 
obtained by requiring that the computed Raa in 0-5% 
central collisions at px = 10 GeV is approximately 0.2. 

Using this q profile, one can compute the medium mod- 
ified fragmentation function D [Eq.Q] for jets which 
originate at r and travel in a given direction 9j(n), decou- 
pling from the medium after traveling a distance Ci(^ 8j) 
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[when T(r + u^lXl) = 160 MeV], with a given initial 
energy p c , and a final hadron momentum specified in 
terms of px and y. The medium modified fragmentation 
function is obtained by solving the differential equation 
in Eq. Q starting with an input vacuum fragmentation 
function at Q 2 — Q 2 min -, and evolving up to Q 2 = p\. 
The physical picture is that of a hard jet that starts with 
a high virtuality Q 2 = p\ and drops down in virtual- 
ly as it splits into less virtual partons while propagating 
through the medium. Eventually at a Q 2 = Q 2 min it exits 
the medium. Thus, there is a strong correlation between 
energy, virtuality and distance travelled by the jet. 

For this attempt to predict the Raa at the LHC, we 
will make two sets of approximations which will include 
the correlation between energy, virtuality and position in 
terms of their averages, thus providing a systematic er- 
ror in extrapolating the calculations at a given px, cen- 
trality and energy of collision to other systems. First, 
we integrate the product of the energy loss kernel K 
in Eq. ([5| and the nuclear overlap function \Taa(5) = 
J d 2 rtA(r — 6/2)ts(r + b/2)] over all jet origins and di- 
rections of exit, for all values of the ratio Q 2 /(2p c ). Note 
that the energy loss kernel correlates location with en- 
ergy p c and virtuality Q 2 through this ratio. We then 
carry out the evolution in virtuality Q 2 for jets with all 
allowed energies using the position and jet direction in- 
tegrated kernel. For each pr of the detected hadron, 
Qmin = (Pc) I ((h) where (p c ) is the mean energy of the 
parent parton in vacuum and ((l) is the mean distance 
travelled by the jets that escape the medium with a virtu- 
ality Q 2 > Q min - For one set of curves (dashed lines in all 
plots) , the mean q of the medium is calculated by averag- 
ing over all possible jet paths in the fluid medium. Then 
a single emission formalism (in a homogeneous static 
medium with the mean q) is used to calculate the length 
(max over which a parton with energy (p c ) will lose all 
its energy. The mean path length ((i) is then calculated 
by averaging over paths in the fluid medium with ( rna x 
as the maximum allowed length. Since a single emission 
is somewhat inefficient in removing the energy from the 
parton, this yields a larger (( L ) leading to smaller Q 2 min 
and thus an over estimate of the quenching. For hadron 
Pt ~ 10 GeV in 0-5% collisions at RHIC, ( max = 5 fm. 
For the other set of curves, (solid lines in all plots) ((l) 
is calculated with the restriction that all lengths larger 
than 5 fm are excluded. Since jets with larger energy 
should penetrate deeper in the medium, this provides an 
underestimate of the quenching. 

With this set up, we can now calculate the invariant 
cross section to produce a high pt hadron in a heavy-ion 
collision using Eq. Q and obtain the Raa using Eq. ([I]) . 
The results for RHIC are shown in Fig. [lj compared with 
the Raa measured by PHENIX in two centrality bins (0- 
5% and 20-30%). As mentioned above, the one dimen- 
sionfull parameter go is dialed so that the calculated Raa 
matches the experimental value in the 0-5% collisions at 



p T = 10 GeV [i.e. R A a(0-5%,Pt = 10GeV) ~ 0.2]. The 
slope of the curve with px and the shift with centrality 
are predictions. The go required is 2.2 GeV 2 /fm. In these 
calculations, Q 2 min = (Pc)/((l) until this value drops be- 
low 1 GeV 2 wherein we hold it fixed at 1 GeV 2 (this 
occurs at a pt ^ 8 GeV). Our choice of go is validated 
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FIG. 1: (Color online) The Raa versus pr for two cen- 
tralities. The green bands represent the calculations pre- 
sented in this Letter with a go — 2.2 GeV 2 /fm. No shad- 
owing is included and the minimum virtuality (Qmin) is se t 
to max{(p c )/((L}, lGeV 2 }. See text for details. 

by comparing with the in-plane and out-of-plane Raa in 
the 20-30% centrality collisions as measured by PHENIX 
(plotted in Fig. [2]) . The solid black curve represents the 
Raa for jets with 6j chosen such that < 3 < 7r/12, for 




a case of a fixed ((„ 



Note, 9j = is the direction of b. 



The solid green line is the Raa with 57r/12 < 9j < 7r/2, 
for a case of a fixed ((max)- The dashed black and 
green lines are the same calculations but with calculated 
((max)- We obtain a good description of the data for 
p T > 8 GeV. For p T < 8 GeV, the Q 2 mm > 1 GeV 2 
restriction begins to affect the calculation. 

With the value of go extracted from RHIC collisions 
and g scaled with the entropy density of the medium 
we now calculate the Raa in Pb-Pb collisions at y/s = 
2.76 ATeV, at the LHC. In order to fit the increased 
charged hadron multiplicities, one requires the fluid dy- 
namical simulation to possess almost twice as much en- 
tropy density at thermalization as at RHIC (T max = 436 
MeV). Also larger velocity gradients drive the system to 
larger sizes than at RHIC prior to freezeout [8j. As a 
result, the mean length traversed by jets escaping the 
medium ((Cl)) is larger. To obtain the Raa we perform 
an almost identical calculation as that for RHIC energies, 
the sole difference being the range of available jet ener- 
gies. As shown in Fig. [3j we obtain a good description of 
the measured Raa (for both 0-5% and 10-30% centrali- 
ties) for pt > 9 GeV. The solid and dashed lines indicate 
the cases of fixed and calculated ((max) as before. This 
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FIG. 2: (Color Online) A parameter free calculation of the 
Raa in plane and out of plane in the 20-30% centrality colli- 
sions and comparison with PHENIX data for the same. 



indicates that, in spite of the produced medium being 
both larger and denser, the basic jet quenching observ- 
ables at high px at the LHC can be described in an iden- 
tical pQCD based formalism as used at RHIC. 
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FIG. 3: (Color Online) A parameter free calculation of the 
Raa in the 0-5% most central and 20-30% central collisions, 
and comparison with CMS data for the same. 

No doubt, future efforts will improve on various ap- 
proximations made in this calculation, in particular the 
treatment of the correlation between energy, virtuality 
and distance travelled by a parton. However it is not clear 
if these corrections are the only source of the discrepancy 
between the calculation and the data at px ~ 6 GeV. 
It is entirely possible that non-perturbative phenomena 
such as recombination extend up to 8-9 GeV at the LHC, 
due to the larger flow at these energies. The extraordi- 
narily high densities may also enhance the magnitude of 
non-linear effects in jet quenching which have so far been 
ignored. In order to gauge this uncertainty, we also plot 



the Raa with a calculated (( ma x) and adjust the qo to 
obtain a better fit with the 0-5% LHC data. This seems 
to require a 30% lower value of go- Another source of 
correction is the drag and diffusion experienced by the 
partons [IB], not included in this calculation 

In conclusion, we have performed a pQCD based cal- 
culation of the suppression of the high px hadron spec- 
tra due to radiative energy loss of hard partons at both 
RHIC and LHC. The local transport coefficient q which 
controls the magnitude of radiative energy loss is scaled 
with the entropy density. The entropy density profile of 
the produced matter is controlled by a 2+1D viscous fluid 
dynamical simulation which has been tuned to describe 
the soft spectrum and elliptic flow. The overall normal- 
ization of q is provided by setting its value at a given T or 
s, in this case qo — 2.2 GeV 2 /fm at the maximum RHIC 
temperature of T = 344MeV. With this value being set, 
we predict the px dependence of the Raa hi 0-5% and 
20-30% centrality events at RHIC as well as the Raa ver- 
sus reaction plane, and the Raa measured in 0-5% and 
20-30% centrality events at ^/s = 2.76 ATeV at the LHC. 
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